PCT 



REQUETE 



Le soussigne requiert que la presente demande 
internationale soit traitee conformement au Traite de 
cooperation en matiere de brevets. 



Reserve a V office recepteur 



Demande internationale n° 



Date du depot international 



Nom de V office recepteur et "Demande internationale PCT" 



Reference du dossier du deposant ou du mandataire (facultatij) 
(12 caracteres au maximum) 4685/00/JC 



Cadre n' I TITRE DE L' INVENTION METHODE POUR MODEL1SER DES DEPLACEMENTS DE 
FLUIDES DANS UN MILIEU POREUX TENANT COMPTE D'EFFETS D'HYSTERESIS 



Cadre n° II DEPOSANT 



Nom et adresse : (Nom de famille suivi du prenom; pour une personne morale, designation 
officielle complete. L 'adresse doit comprendre le code postal et le nom du pays. Le pays de 
I adresse indiquee dans ce cadre est I 'Etat ou le deposant a son domicile si ducun domicile 
n 'est indique ci-dessous.) 

INST1TUT FRANCAIS DU PETROLE 
1 & 4 avenue de Bois Preau 
92852 Rueil-Malmaison cedex 
FRANCE 



□ 



Cette personne est aussi 
inventeur. 



n° de telephone 

01 47 52 60 00 



n° de teiecopieur 

01 47 52 70 03 



n° de teleimprimeur 



Nationality (nom de TEtat) : 
FRANCE 



Domicile (nom de l'Etat) 
FRANCE 



Cette personne est 
deposant pour : 



□ tous les Etats 
designes 



[TTl tous, les Etats designes sauf I I les Etats-Unis d'Amerique X~\ les Etats indiques dans 
1 A I les Etats- Unis d'Amerique | | seulement | | l e cadre supplementaire 



Cadre n° III AUTRE(S) DEPOSANT(S) OU (AUTRE(S)) INVENTEUR(S) 



Nom et adresse : (Nom de famille suivi du prenom; pour une personne morale, designation 
officielle complete. L 'adresse doit comprendre le code postal et le nom du pays. Le pays de 
radresse indiquee dans ce cadre est I 'Etat ou le deposant a son domicile si aucun domicile 
n 'est indique ci-dessous.) 

GAZ DE FRANCE 

23 rue Philibert Delorme 

75840 Paris cedex 1 7 

FRANCE 


Cette personne est : 

| y | deposant seulement 

[ | deposant et inventeur 

j | inventeur seulement 
(Si cette case est cochee, 
ne pas remplir la suite.) 


National ite (nom de PEtat) : 
FRANCE 


Domicile (nom de FEtat) : 
FRANCE 



Cette personne est 
deposant pour : 



I | tous les Etats rZjr] tous les Etats designes sauf I I les Etats-Unis d'Amerique I ] les Etats indigues 
I I designes I /» I les Etats-Unis d'Amerique | _| seulement | | le cadre supplemer 



dans 
iementaire 



| Xl D'autres deposants ou inventeurs sont indiques sur une feuille annexe. 



Cadre n° IV MANDATAIRE OU REPRESENTANT COMMUN; OU ADRESSE POUR LA CORRESPONDANCE 



La personne dont Pidentite est donnee ci-dessous est/a ete designee pour agir au nom du ou | 1 mo „ a~*~~ rm 

des deposants aupres des autorites intemationales com P aentes%oinme: ■ L_| mandataire LXJ representant commun 



Nom et adresse : (Nom de famille suivi du prenom; pour une personne morale, desigt 
1 complete. L 'adresse doit comprendre le code postal et le nom du pays.) 

INSTITUT FRANCAIS DU PETROLE 


lation officielle 


n° de telephone 
01 47 52 60 00 


1 & 4 avenue de Bois Preau 
92852 Rueil-Malmaison cedex 
FRANCE 




n° de teiecopieur 

01 47 52 70 03 




n° de teleimprimeur 



□ Adresse pour la correspondance : cocher cette case lorsque aucun mandataire ni representant commun n'est/n'a ete designe 
et que Tespace ci-dessus est utilise pour indiquer une adresse speciale a laquelle la correspondance doit etre envoyee. 



Formulaire PCT/RO/101 (premiere feuille) Quillet 1998; reimpression janvier 2001) 



Voir les notes relatives au formulaire de requite 



m 
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Feuille n° . 



Suite du cadre n° III AUTRE(S) DEPOSANT(S) OU (AUTRE(S)) INVENTEUR(S) 


Si aucun des sous-cadres suivants n 'est utilise, cette feuille ne doit pas itre incluse dans la requite. 


Nom et adresse : (Nom de famille suivi du prenom; pour une personne morale, designation 
officielle complete. L 'adresse doit comprendre le code postal et le nom du pavs. Le pavs de 
F adresse indiquee dans ce cadre est I Etat oil le deposant a son domicile si aucun domicile 
n est indique ci-dessous.) 

EGERMANN Patrick 
193 route de I'Empereur 
92500 Rueil-Malmaison 
FRANCE 


Cette personne est : 

| | deposant seuiement 

|y | deposant et inventeur 

[ "J inventeur seuiement 
(Si cette case est cochee, 
ne pas remplir la suite.) 


Nationality (nom de TEtat) : 
FRANCE 


Domicile (nom de TEtat) : 
FRANCE 


Cette personne est I 1 tousles Etats I 1 tousles Etats designes sauf les Etats-Unis d'Amerique 1 I les Etats indiques dans 
deposant pour : 1 1 designes 1 1 les Etats-Unis d'Amerique U5u seuiement 1 1 le cadre suppiementaire 


Nom et adresse : (Nom de famille suivi du prenom; pour une personne morale, designation 
officielle complete. L 'adresse doit comprendre le code postal et le nom du pays. Le pays de 
radresse indiquee dans ce cadre est I 'Etat oil le deposant a son domicile si aucun domicile 
n 'est indique ci-dessous.) 

VIZIKA-KAVVADIAS Olga 

3, allee de la Pagerie 

92500 Rueil-Malmaison 

FRANCE 


Cette personne est : 

| | deposant seuiement 

|X | deposant et inventeur 

| . | inventeur seuiement 

(Si cette case est cochee, • 
ne pas remplir la suite.) 


Nationality (nom de TEtat) : 
FRANCE 


Domicile (nom de TEtat) : 
FRANCE 


Cette personne est 1 1 tousles Etats 1 | tous.les Etats designes sauf pZT] les Etats-Unis d'Amerique 1 1 les Etats indiques dans 

deposant pour : 1 1 designes | | les Etats-Unis d'Amerique 1 A 1 seuiement | | | e cadre suppiementaire 


Nom et adresse : (Nom de famille suivi du prenom; pour une personne morale, designation 
officielle complete. L 'adresse doit comprendre le code postal et le nom du pays. Le pays de 
radresse indiquee dans ce cadre est I 'Etat oil le deposant a son domicile si aucun domicile 
n 'est indique ci-dessous.) 

DALLET Laurent 
57 quai de Grenelle 
75015 Paris 
FRANCE 


Cette personne est : 

| J deposant seuiement 

|y | deposant et inventeur 

| | inventeur seuiement 
(Si cette case est cochee. 
ne pas remplir la suite.) 


Nationality (nom de TEtat) : 
FRANCE 


Domicile (nom de TEtat) : | 
FRANCE | 


Cette personne est I | tous les Etats [ 1 tous.les Etats designes sauf ^ les Etats-Unis d'Amerique 1 1 les Etats indiques dans 

deposant pour : 1 1 designes | J les Etats-Unis d'Amerique A seuiement | | le cadre suppiementaire 


Nom et adresse : (Nom de famille suivi du prenom; pour une personne morale, designation 
officielle complete. L 'adresse doit comprendre le code postal et le nom du pays. Le pays de 
radresse indiquee dans ce cadre est I 'Etat oit le deposant a son domicile si aucun domicile 
n 'est indique ci-dessous.) 

KALAYDJ1AN Frangois 

1 7 rue Danton 

92500 Rueil-Malmaison 

FRANCE 


Cette personne est : 

| | deposant seuiement 

|X| deposant et inventeur I 

[^] inventeur seuiement 
(Si cette case est cochee, 
ne pas remplir la suite.) 


Nationality (nom de TEtat) : 
FRANCE 


Domicile (nom de TEtat) : 
FRANCE 


Cette personne est | I tous les Etats | I tous.les Etats designes sauf r^l les Etats-Unis d'Amerique 1 1 les Etats indiques dans 

deposant pour : 1 1 designes | | les Etats-Unis d'Amerique IA I seuiement | | l e cadre suppiementaire 


ly| D'autres deposants ou inventeurs sont indiques sur une autre feuille annexe. 



Formulaire PCT/RO/10 1 (feuille annexe) (juillet 1 998; reimpression janvier 200 1 ) Voir les notes relatives au formulaire de requite 
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Suite du cadre n° III AUTRE(S) DEPOSANT(S) OU (AUTRE(S)) INVENTEUR(S) 


Si aucun des sous-cadres suivants n'est utilise, cette feuille ne doit pas etre bid use dans la requite* 


Nom et adresse : (Norn de famille snivi du prenom; pour une personne morale, designation 
officielle complete. L 'adresse doit comprendre le code postal et le nom du pays. Le pays de 
F adresse indiquee dans ce cadre est I 'Etat oit le deposant a son domicile si ducun domicile 
n est indique ci-dessous.) 

REQUIN Christophe 

1 9 rue du Perou 

76000 Rouen 

FRANCE 


Cette personne est : j 
f" [ deposant seulement 

\}C\ deposant et inventeur 

[ ~] inventeur seulement 
(Si cette case est cochee, 
ne pas remplir la suite.) 


Nationalite (nom de PEtat) : 
FRANCE 


Domicile (nom de PEtat) : 
FRANCE 


Cette personne est |— | tousles Etats | 1 tousles Etats designes sauf y, les Etats-Unis d'Amerique 1 1 les Etats indiques dans 
deposant pour : 1 1 designes 1 1 les Etats-Unis d'Amerique l*J seulement 1 I le cadre supplemental 


Nom et adresse : (Nom de famille suivi du prenom; pour une personne morale, designation 
officielle complete. L 'adresse doit comprendre le code postal et le nom du pays, Le pavs de 
I adresse indiquee dans ce cadre est VEtat ou le deposant a son domicile si aucun domicile 
n 'est indique ci-dessous.) 


Cette personne est : 
\~\ deposant seulement 

j ~| deposant et inventeur 

[ J inventeur seulement 
(Si cette case est cochee, 
ne pas remplir la suite.) \ 


Nationality (nom de I 'Etat) : 


Domicile (nom de PEtat) : 


Cette personne est | 1 tousles Etats I | tousjes Etats designes sauf 1 | les Etats-Unis d'Amerique | j les Etats indiques dans 

deposant pour : 1 1 designes | J ies Etats-Unis d'Amerique | | seulement | | le cadre supplementaire 


Nom et adresse : (Nom de famille suivi du prenom; pour une personne morale, designation 
officielle complete. L 'adresse doit comprendre le code postal et le nom du pays. Le pays de 
V adresse indiquee dans ce cadre est I 'Etat oit le deposant a son domicile si aucun domicile 
n est indique ci-dessous.) 


Cette personne est : 
[~] deposant seulement 

deposant et inventeur 

j | inventeur seulement 
(Si cette case est cochee. 
ne pas remplir la suite.) 


Nationalite (nom de PEtat) : 


Domicile (nom de PEtat) : 


Cette personne est 1 1 tous les Etats I I tous. les Etats designes sauf 1 1 les Etats-Unis d'Amerique I 1 les Etats indiques dans 

deposant pour : 1 1 designes | | | e s Etats-Unis d'Amerique | | seulement J | le cadre suppfementaire 


Nom et adresse : (Nom de famille suivi du prenom; pour une personne morale, designation 
officielle complete. L adresse doit comprendre le code postal et le nom du pays. Le pays de 
I adresse indiquee dans ce cadre est VEtat ou le deposant a son domicile si aucun domicile 
n est indique ci-dessous.) 


Cette personne est : 

| | deposant seulement 

[~] deposant et inventeur 

| [ inventeur seulement 
(Si cette case est cochee, 
ne pas remplir la suite.) 


Nationalite (nom de PEtat) : 


Domicile (nom de PEtat) : 


Cette personne est | 1 tous les Etats I I tous, les Etats designes sauf 1 | les Etats-Unis d'Amerique 1 1 les Etats indiques dans 

deposant pour : 1 I designes | | les Etats-Unis d'Amerique | | seulement | | l e cadre supplementaire 


1"""*] D'autres deposants ou inventeurs sont indiques sur une autre feuille annexe. 



Formulaire PCT/RO/1 0 1 (feuille annexe) (juillet 1 998; reimpression janvier 200 1 ) Voir les notes relatives auformulaire de requete 
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Cadre n° V DESIGNATION D'E 




Feuille if 




• Les designations suivantes sont faites conformement a la regie 4. 9. a) (cocher les cases appropriees; une au moins doit I 'etre) : 
Brevet regional 

Q AP Brevet ARIPO : GH Ghana, GM Gambie, KE Kenya, LS Lesotho, MVV Malawi, MZ Mozambique, SD Soudan, 

SL Sierra Leone, SZ Swaziland, TZ Republique-Unie de Tanzanie, UG Ouganda, Z\V Zimbabwe et tout autre Etat qui est un Etat 

contractant du Protocole de Harare et du PCT 
0 EA Brevet eurasien : AM Armenie, AZ Azerbaidjan, BY Belarus, KG Kirghizistan, KZ Kazakhstan, MD Republique de Moldova. 

RU Federation de Russie, TJ Tadjikistan. TM Turkmenistan et tout autre Etat qui est un Etat contractant de la Convention sur 

le brevet eurasien et du PCT 

0 EP Brevet europeen : AT Autriche, BE Belgique, CH et LI Suisse et Liechtenstein, CY Chypre, DE Allemasne. 
DK Danemark, ES Espagne, FI Finlande, FR France, GB Royaume-Uni, GR Grece, IE Irlande, IT Italie, 
LU Luxembourg, MC Monaco, NL Pays-Bas, PT Portugal, SE Suede, TR Turquie et tout autre Etat qui est un Etat contractant 
de la Convention sur le brevet europeen et du PCT 

0 OA Brevet OAPI : BF Burkina Faso, BJ Benin, CF Republique centrafricaine, CG Congo, CI Cote d'lvoire, 
CM Cameroun, GA Gabon, GN CJuinee, GW Guinee-Bissau, ML Mali, MR Mauritanie, NE Niger, SN Senegal, 
TD Tchad, TG Togo et tout autre Etat qui est un Etat membre de l'OAPI et un Etat contractant du PCT (si une autre forme 
de protection ou de traitement est souhaitee, le preciser sur la ligne pointillee) 

Brevet national (si une autre forme de protection ou de traitement est souhaitee, le preciser sur la ligne pointillee) : 

0 AE Emirats arabes unis 0 LC Sainte-Lucie 

0 AG Antigua-et-Barbuda 0 LK Sri Lanka 

0 AL Albanie 0 LR Liberia 

STI AM Armenie 0 LS Lesotho 

0 AT Autriche 0 LT Lituanie 

0 AU Austral ie 0 LU Luxembourg 

0 AZ Azerbaidjan ' 0 LV Lettonie 

0 BA Bosnie-Herzegovine 0 MA Maroc 

0 BB Barbade 0 MD Republique de Moldova 

0 BG Bulgarie 0 MG Madagascar 

0 BR Bresil 0 MK Ex-Republique yougoslave de Macedoine 

0 BY Belarus 0 MN Mongolie 

□ BZ Belize 0 MW Malawi 

0 CA Canada 0 MX Mexique 

0 CH et LI Suisse et Liechtenstein 0 MZ Mozambique 

0 CN Chine . . 0 NO Norvege 

0 CR Costa Rica □ NZ. Nouvelle-Zelande 

0 CU Cuba 0 PL Pologne 

0 CZ Republique tcheque 0 PT Portugal 

0 DE Allemagne 0 RO Roumanie 

0 DK Danemark . . . . 0 RU Federation de Russie 

0 DM Dominique 0 SD Soudan 

0 DZ Algerie 0 SE Suede 

0 EE Estonie Q SG Singapour 

0 ES Espagne 0 SI Slovenie 

0 FI Finlande 0 SK Slovaquie 

0 GB Royaume-Uni 0 SL Sierra Leone 

0 GD Grenade 0 TJ Tadjikistan 

0 GE Georgie 0 TM Turkmenistan 

0 GH Ghana 0 TR Turquie 

0 GM Gambie 0 TT Trinite-et-Tobago 

0 HR Croatie 0 TZ Republique-Unie de Tanzanie 

0 HU Hongrie 0 UA Ukraine 

0 ID Indonesie 0 UG Ouganda 

0 IL Israel 0 US Etats-Unis d'Amerique 

0 IN Inde □ UZ Ouzbekistan 

0 IS Islande 0 VN Viet Nam 

0 JP Japon 0 YU Yougoslavie 

0 KE Kenya 0 ZA Afrique du Sud 

0 KG Kirghizistan 0 ZW Zimbabwe 



Declaration concernant les designations de precaution : outre les designations faites ci-dessus, le deposant fait aussi conformement 
a la regie 4.9.b) toutes les designations qui seraient autorisees en vertu du PCT, a 1 'exception de toute designation indiquee dans le cadre 
supplemental comme etant exclue de la portee de cette declaration. Le deposant declare que ces designations additionnelles sont 
faites sous reserve de confirmation et que toute designation qui n'est pas confirmee avant V expiration d'un delai de 1 5 mois a compter 
de la date de priorite doit etre consideree comme retiree par le deposant a l'expiration de ce de"lai. (La confirmation (y compris les taxes) 
doit parvenir a I 'office recepteurdans le delai de 1 5 mois.) 



0 KP 
□ KR 
0 KZ 



Republique populaire democratique de Coree 

Republique de Coree 

Kazakhstan 



Case reservee pour la designation d' Etats qui sont devenus parties au 
PCT apres la publication de la presente feuille : 

0 



Formulaire PCT/RO/101 (deuxieme feuille) (Janvier 2001) 



Voir les notes relatives au formulaire de requite 
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Cadre n" VI 



Feuille n° 



REVENDICATION DE PRIORITE 



□ 



T^^tres 



D autres revendications de priorite sont 
indiquees dans le cadre supplementaire. 



Date de depot 
de la demande anterieure 
(jour/mois/annee) 


Numero 
de la demande anterieure 


Lorsque la demande anterieure est une : 


demande nationale : 
pays 


demande regionale :* 
office regional 


demande Internationale : 
office recepteur 


0) 17juillet2000 
(17/07/00) 


00/09368 


FRANCE 






(2) 










(3) 











| x I L' office recepteur est prie de preparer et de transmettre au Bureau international une copie certifiee conforme de la ou des demandes 
anterieures (settlement si la demande anterieure a ete deposee aupres de I 'office qui, aux Jins de nn/nQ Q C Q 
la presente demande Internationale, est I 'office recepteur) indiquees ci-dessus au(x) point(s) : 00/09368 



* Si la demande anterieure est une demande ARIPO, il est obligatoire d'indiquer dans le cadre supplementaire au moins un nays partie a la Convention 
de Paris pour la protection de la propriete industrielle pour lequel cette demande anterieure a ete deposee (regie 4.10.b)ii)). Voir le cadre supplementaire. 



Cadre n w VII ADMINISTRATION CHARGEE DE LA RECHERCHE INTERNATIONALE 



Choix de Fadministration chargee de la recherche 
internationale (ISA) (si plusieurs administrations 
chargees de la recherche internationale sont competentes 
pour proceder a la recherche internationale, indiquer 
1' administration choisie; le code a deux lettres peut etre 
utilise) : 

ISA/ 



Demande d'utilisation des resultats d'une recherche anterieure; mention de 
cette recherche (si une recherche anterieure a ete effectuee par I 'administration 
chargee de la recherche internationale ou demandee a cette dernier e) : 

Date (jour/mois/annee) Numero Pays (ou office regional) 



13 juin 2001 



FA 592697 



FRANCE 



Cadre n tt VIII BORDEREAU; LANGUE DE DEPOT 



La presente demande internationale contient 
le nombre de feuilles suivant : 



requete 

description (sauf partie reservee 
au listage des sequences) 

revendications 

abrege 

dessins 



18 
2 
1 
5 



partie de la description reservee 
au listage des sequences : 

Nombre total de feuilles 



"31 



Le ou Ies elements coches ci-apres sont joints a la presente demande internationale : 

1 . H feuille de calcul des taxes 

2. □ pouvoir distinct signe 

3. □ copie du pouvoir general; numero de reference, le cas echeant : 

4. □ explication de V absence d'une signature 

5. □ document(s) de priorite indique(s) dans le cadre n° VI au(x) point(s) : 

6. □ traduction de la demande internationale en (langue) : 

7. □ indications separees concernant des micro-organismes ou autre materiel 

biologique deposes 

8. □ listage des sequences de nucleotides ou d'acides amines sous forme 

dechiffrable par ordinateur 

9. □ autres elements (precis er) : copie rapport de recherche franca is 



Figure des dessins qui ~* 
doit accompagner V abrege : ^ 


Langue de depot de la ! 
demande internationale : frangaiS 


Cadre n° IX SIGNATURE DU DEPOSANT OU DU MANDATAIRE 



A cote de chaque signature, indiquer le nom du signataire et, si cela n 'apparait pas clairement a la lecture de la requete, a quel litre I inter esse signe. 

Pour Alfred ELMALEH, 



^feahCOADOUR 
(jngenieur 



Reserve a l'office recepteur 



1 . Date effective de reception des pieces supposees 
constituer la demande internationale : 



3. Date effective de reception, rectifiee en raison de la reception ulte- 
rieure, mais dans les delais, de documents ou de dessins completant ce 
• qui est suppose constituer la demande internationale : 



4. Date de reception, dans les delais, des corrections 
demandees selon 1' article 1 1.2) du PCT : 



2. Dessins : 
| | recus : 



□ 



non recus : 



5. Administration chargee de la recherche 
internationale (si plusieurs sont competentes) : 



ISA/ 



□ Transmission de la copie de recherche differee 
jusqu'au paiement de la taxe de recherche. 



Reserve au Bureau international 



Date de reception de l'exemplaire 
original par le Bureau international : 



Formulaire PCT/RO/101 (derniere feuille) Quillet 1998; reimpressionjanvier 2001) 



Voir les notes relatives au formulaire de requete 



41 < • 



This Page Blank (uspto) 



Cette feuille ne fait pas po^^^m la demande Internationale ni ne compte con^^Bie feuille de celle-ci 



PCT 



FEUILLE DE CALCUL DES TAXES 
Annexe de la requete 



Reference du dossier du 

deposant ou du mandataire 4685/00/JC 



Reserve a 1'office recepteur 



Demande internaiionale n u 



Timbre a date de 1'office recepteur 



Deposant 



INSTITUT FRANCAIS DU PETROLE et GAZ DE FRANCE 



CALCUL DES TAXES PRESCRITES 



TAXE DE TRANSMISSION 



400,00 



2. TAXE DE RECHERCHE 

Recherche internaiionale a effectuer par 

(Si plusieurs administrations chargees de la recherche Internationale sont 
competentes en ce qui concerne la demande Internationale, inscrire le nom de celle 
qui est choisie pour la recherche inter nationale.) 

3. TAXE INTERNATIONALE 



6198,79 



Taxe de base 

La demande Internationale contient 31 



59,04 



30 premieres feuilles 
1 

feuilles suivantes 

Additionner les montants portes dans les cadres 
bl et b2 et inscrire le total dans le cadre B . . . 



feuilles. 

I 2682,66 



bi 



montant additionnel 



= | 59,04 



b2 



Taxes de designation 

La demande intemationale contient 87 



designations. 



2741.90 



577.24 



nombre de taxes de 
designation dues (maximum 6) 



montant de la taxe de designation 



= | 3463. 



44 



Additionner les montants portes dans les cadres B et D, et 

inscrire le total dans le cadre I 

(Les deposants de certains Etats ont droit a une reduction de 75 % sur la taxe 
Internationale. Lorsque le deposant a (ou tous les deposants ont) droit a cette 
reduction, la somme devant Jigurer sous 1 est egale a 25 % de la somme des 
montants figurant sous B et D.) 

TAXE AFFERENTE AU DOCUMENT DE PRIORITE (le cas echeant) . 

TOTAL DES TAXES DUES 

Additionner les montants portes dans les cadres 

T, S, I et P, et inscrire le resultat dans le cadre TOTAL 



6205.34 



100 



12904.13 



TOTAL 



! ] Les taxes de designation seront payees ulterieurement. 



MODE DE PAIEMENT 

I x I autorisation de debiter un compte l~~ | traite bancaire 

L— > de depot (voir ci-dessous) ^—^ 

□ 

cheque Q 

\ ] mandat postal |" | 



especes 

timbres fiscaux 



coupons 
| | autres (preciser): 



AUTORISATION CONCERNANT UN COMPTE DE DEPOT (les offices recepteurs ne permettent pas tous lutilisation de ce mode de paiement) 



L* office recepteur/ 



□ 
□ 



239 



Numero du compte de depot 



est autorise a debiter mon compte de depot du total des taxes indique ci-dessus. 

(cette case ne peut etre cochee que si les conditions relatives aux comptes de depot etablies par I 'office 
recepteur le permettent) est autorise a debiter mon compte de depot de tout montant manquant - ou a le 
crediter de tout excedent - dans le paiement du total des taxes indique ci-dessus. 

est autorise a debiter mon compte de depot du montant de la taxe afferente a l'etablissement du document 
de priorite et a sa transmission au Bureau international de I'OMPI. 



09/07/00 



Date (jour/mois/annee) 



Signature 



Jean COADOUR 



Formulaire PCT/RO/101 (Annexe) (janvier 2001) 



Voir les notes relatives a la feuille de calcul des taxes 



\ 
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TRAITE DE C 



ERATION EN MATIERE DE B 

PCT 



RAPPORT DE RECHERCHE INTERNATIONALE 
(article 18 et regies 43 et 44 du PCT) 



Reference du dossier du deposant ou 
du mandataire 

4685/00/JC 


POUR SUITE volr la notification de transmission du rapport de recherche internationale 
(formulaire PCT/IS A/220) et, le cas echeant, le point 5 ci-apres 

A DONNER 


Demande Internationale n° 

PCT/FR 01/02212 


Date du depot international (jour/mois/annee) 

10/07/2001 


(Date de priorite (la plus ancienne) 
(jour/mois/annee) 

17/07/2000 


Deposant 

INSTITUT FRANCAIS DU PETROLE 



Le present rapport de recherche internationale, etabli par I'administration chargee de la recherche internationale, est transmis au 
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1 . Base du rapport 

a. En ce qui concerne la langue, la recherche internationale a ete effectuee sur la base de la demande internationale dans la 
langue dans laquelle elle a ete deposee, sauf indication contraire donnee sous le meme point. 



□ 



la recherche internationale a ete effectuee sur la base d'une traduction de la demande internationale remise a I'administration. 



En ce qui concerne les sequences de nucleotides ou d'acides amines divulguees dans la demande internationale (te cas echeant), 
la recherche internationale a ete effectuee sur la base du I i stage des sequences : 

[ | contenu dans la demande internationale, sous forme ecrite. 

deposee avec la demande internationale, sous forme dechiffrable par ordinateur. 
remis ulterieurement a I'administration, sous forme ecrite. 
remis ulterieurement a I'administration, sous forme dechiffrable par ordinateur. 



□ 
□ 
□ 
□ 

□ 



La declaration, seion laquelie le listage des sequences presente par ecrit et fourni ulterieurement ne vas pas au-dela de ta 
divulgation faite dans la demande telle que deposee, a ete fournie. 

La declaration, selon laquelle les informations enregistrees sous forme dechiffrable par ordinateur sont identiques a celles 
du listage des sequences presente par ecrit, a ete fournie. 



2. 
3. 



| | II a ete estime que certaines revendications ne pouvaient pas faire I'objet d'une recherche (voir le cadre I). 
I ] II y a absence d 'unite de I'invention (voir le cadre II). 



4. En ce qui concerne le titre, 

[X] le texte est approuve tel qu'il a ete remis par le deposant. 

[ | Le texte a ete etabli par I'administration et a la teneur suivante: 



5. En ce qui concerne I'abrege, 

j | le texte est approuve tel qu'il a ete remis par le deposant 

ry-i le texte (reproduit dans le cadre III) a ete etabli par I'administration conformement a la regie 38.2b). Le deposant peut 
LAJ presenter des observations a I'administration dans un delai d'un mois a compter de la date cfexpedition du present rapport 
de recherche internationale. 

6. La figure des dessins a publier avec I'abrege est la Figure n° 3 

[X] suggeree par le deposant. [_J Aucune des figures 

□~ , , 4 , , r n'est a publier. 

parce que le deposant n'a pas suggere de figure. 

[ J parce que cette figure caracterise mieux I'invention. 
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Cadre III TEXTE DE L'ABREGE (suite du point 5 de la premiere feuille) 



La methode comporte la determination experimentale de la courbe de variation 

de la pression capillaire dans les pores en fonction de la saturation 

relativement aux phases liquides, la modelisation des pores du milieu poreux 

par une distribution de capillaires a repartition fractale en considerant, 

dans le cas d'un melange triphasique eau (fluide mouillant) -huile-gaz par 

exemple une stratification des constituants a l'interieur des pores, 

avec l'eau au contact des parois, le gaz au centre et l'huile formant 

une couche intercalaire , la determination, a partir de cette courbe de 

pression capillaire, des valeurs de dimension fractale correspondant a 

une serie de valeurs donnees de la saturation relativement a la phase 

liquide, la modelisation des effets d' hysteresis modifiant les 

saturations mobiles des fluides ef f ectivement deplaces dans 

1 ' echantillon qui varient au cours de cycles de drainage et 

d' imbibition, la modelisation des permeabilites relatives directement 

sous forme d' expressions analytiques dependant des dif f erentes valeurs 

de dimension fractale obtenues et en accord avec la repartition 

stratifiee des differents fluides dans les pores, et 1 ' utilisation d'un 

simulateur de milieu poreux pour determiner, a partir des permeabilites 

relatives, les conditions optimales de Replacement des fluides dans le 

mi 1 ieu poreux . 
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(57) Abstract: The invention concerns a method comprising experimentally deter- 
mining the variation curve of capillary pressure in pores on the basis of saturation 
with respect to liquid phases, modelling the pores of the porous environment by 
capillary distribution with fractal distribution while taking into account in the case 
of a three-phase water (wetting fluid)-oil-gas mixture for example stratification of 
the constituents inside the pores, with water in contact with the walls, gas in the 
centre and oil forming an intervening Jayer, determining, on the basis of said capil- 
lary pressure curve, values of fractal dimension corresponding to a series of given 
values of saturation relative to the liquid phase, modelling hysteresis effects modi- 
fying the mobile saturation levels of fluids actually displaced in the sample which 
vary in the course of drainage and imbibition cycles, modelling relative permeabil- 
ity levels directly in the form of analytical expressions depending on the different 
values of fractal dimension obtained and in accordance with the stratified distribu- 
tion of the different fluids in the pores, and using a porous medium simulator to 
determine, on the basis of the relative permeability levels, the optimal conditions 
of fluid displacement in a porous medium. 



W... WETTING FLUID 
G . GAS PHASE 
RG. .RADIUS (PORE VOLUME) 
O... INTERVENING LAYER 



(57) Abrege : La methode comporte la determination experimentale de la courbe 
de variation de la pression capillaire dans les pores en fonction de la saturation re- 
lativement aux phases liquides, la modelisation des pores du milieu poreux par une 
distribution de capillaires a repartition fractale en considerant, dans le cas d'un me- 
lange triphasique eau (fluide mouillant)-huile-gaz par exemple une stratification des constituants a Tinterieur des pores, avec 1'eau 
au contact des parois, le gaz au centre et 1'huile formant une couche intercalate, la determination, a partir de cette courbe de pression 
capillaire, des valeurs de dimension fractale correspondant a une serie de valeurs donnees de la saturation relativement a la phase 
liquide, la modelisation des effets d'hysteresis modifiant les saturations mobiles des fluides effectivement deplaces dans l'echan- 
tillon qui varient au cours de cycles de drainage et d 7 imbibition, la modelisation des permeabilites relatives directement sous forme 
depressions analytiques dependant des diffcrentes valeurs de dimension fractale obtenues et en accord avec la repartition stratifiee 
des differents fluides dans les pores, et V utilisation d'un simulaleur de milieu poreux pour determiner, a partir des permeabilites 
relatives, les conditions optimales de deplacement des fluides dans le milieu poreux. 
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5 METHODE POUR MODELISER DES DEPLACEMENTS DE FLUIDES 
DANS UN MILIEU POREUX TENANT COMPTE D'EFFETS 
D'HYSTERESIS 

La presente invention concerne une methode pour modeliser des ecoulements 
di ou triphasiques dans un milieu poreux, en drainage et en imbibition. Elle est basee 
10 sur une representation fractale du milieu poreux et sur une approche originale pour 
traiter les phenomenes lies a Thysteresis (changement du sens de variation des 
saturations). 



Ton veut modeliser les ecoulements de fluides dans des milieux poreux, dans le but 
15 d'optimiser les conditions de leur deplacement en drainage et en imbibition. Parmi les 
domaines d'application, on peut citer par exemple : 

a) l'exploitation d'un gisement petrolier et notamment la production assistee 
d'hydrocarbures par injection de fluides, en utilisant par exemple des injections 
alternees de bouchons de liquide et de gaz (methode dite WAG). Elle constitue 
20 pour l'ingenieur de gisements, un bon outil pour etudier egalement les problemes 

de productivite et d'injectivite des puits. 



La methode selon Tinvention est applicable dans de nombreux domaines ou 



b) b) la depollution des sols et notamment des sites industriels, par injection dans 
des couches polluees, de substances telles que des tensio-actifs ; 
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c) le nettoyage de cakes de filtration par deplacement des substances qui y sont 
retenus ; 

d) le sechage du bois ; 

e) Toptimisation de reactions chimiques par deplacement par exemple de produits 
5 reactionnels dans une masse de catalyseur, pour augmenter les surfaces de 

contact, etc. 

La methode selon 1'invention est directement applicable par les ingenieurs de 
reservoir *p<3ur determiner par exemple la methode de recuperation assistee la plus 
appropriee a appliquer a un gisement souterrain d ! hydrocarbures. La methode peut 
10 egalement servir dans le cadre d'operations de depollution de terrains industriels par 
exemple. 

Etat de la technique 

1) Etudes experimentales : 

La determination experimentale des permeabilites relatives d'un milieu 
15 poreux ou se deplace un fluide polyphasique n'est pas une tache aisee. On simplifie 
habituellement les operations de mesure en considerant que Tune des phases est 
immobile dans un etat de saturation irreductible. 

Les valeurs sont acquises par exemple par une mdthode experimentale connue 
dite de "steady state" pour determiner les permeabilites relatives et qui consiste h. 

20 laisser s f ecouler un fluide triphasique avec des rapports de debit imposes entre les 
phases. Les permeabilites relatives exprimees en fonction des deux saturations, sont 
calculees en appliquant la loi de Darcy a chaque phase. H est loin d'etre etabli que les 
mesures des permeabilites relatives obtenues par cette methode soient bien 
representatives des deplacements des fluides et, de toute fagon, elles sont longues a 

25 obtenir car a chaque changement de regime, il faut attendre qu ! un regime d f equilibre 
s'etablisse. 
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Une autre methode connue consiste h realiser des essais en laboratoire pour 
determiner des tableaux de mesures (tel que celui de la Fig.l) reliant les 
permeabilites relatives et les saturations pour chaque couple de fluides du melange 
triphasique. Par ajustement des courbes de production experimentales, on essaie 
5 progressivement d'ajuster les permeabilites relatives triphasiques. Ces tableaux de 
donnees sont introduits alors dans un simulateur numerique de type Athos® qui 
calcule les productions de fluides. Cette methode etant basee sur l'acquisition 
prealable de nombreuses mesures experimentales ajustees progressivement par 
calage, est longue a mettre en ceuvre. 

10 2) Modeles de permeabilites relatives 

Le modele empirique connu dit de Stone permet, par des correlations 
empiriques, de predire des donnees relatives a un ecoulement triphasique a partir de 
donnees correspondant a un ecoulement diphasique. H est valide seulement dans les 
cas d'une forte mouillabilite h Teau et on le considere generalement comme 
15 f aiblement predictif . 

On connait deux types de modeles physiques pour modeliser les ecoulements 
triphasiques, bases sur les courbes de pression capillaire. Les courbes de pression 
capillaire sont en relation avec une saturation (par exemple celle du mercure injecte) 
et un rayon de pore, auquel s f arrete le mercure pour une pression d'injection donnee, 

20 determinee par la loi de Laplace, Pinj.= — 

r 

Un premier modele de representation des milieux poreux, est decrit par : 

- Burdine, N.T.: "Relative permeability calculations from pore size distribution 
data" Trans AIME (1953) Vol. 198 ; ou par 

- Corey, A.T.: "The interrelation between oil and gas relative permeabilities", Prod 
25 Monthly (1954) Vol. 19,38. 
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Suivant ce modele, le milieu poreux est represents par un faisceau de 
capillaires cylindriques avec une distribution de rayons qui sont donnes par la courbe 
de pression capillaire obtenue par injection de mercure. Les permeabilites sont 
obtenues par application de la loi de Poiseuille a Tecoulement de fluides dans ces 
5 capillaires. 

Ce modele est base sur la representation du milieu poreux comme un 
assemblage de capillaires de rayons differents. La relation entre le volume et le rayon 
des pores est donnee par la valeur de la pente du pseudo-plateau. Les trois fluides 
sont censes se partager les capillaires, le fluide mouillant (Teau) occupant les plus 

10 petits, le fluide le moins mouillant (le gaz), les plus gros, le troisieme fluide (Thuile) 
une zone de pores de dimensions intermediaries, n n'est pas possible de decrire les 
interactions entre les fluides puisque dans un tel modele, ils s'ecoulent dans des 
chenaux separes. Ce modele enfin, ne peut etre utile que si le pseudo-plateau couvre 
une large gamme de saturations. Suivant ce modele, les trois phases d'un ecoulement 

15 triphasique se deplacent dans des capillaires differents et il n f y a entre elles aucune 
interaction. 

Un autre mod&le physique connu de representation d ! un milieu poreux est 
decrit par : 

- de Gennes, P.G.: "Partial Filling of a Fractal Structure by a Wetting Fluid" 
20 Physics of Disordered Materials 227-241, New York Plenum Pub. Corp. (1985) ; 

repris par 

- Lenormand, R. "Gravity Assisted Inert Gas Injection: Micromodel Experiments 
and Model based on Fractal Roughness", The European Oil and Gas Conference 
Altavilla Milica, Palermo, Sicily, (1990). 

25 Suivant ce modele, on considere que la surface interne des pores est isotrope 

et a caractere fractal et peut etre modelisee sous la forme d'un "bouquet" de rainures 
capillaires paralleles qui font que les pores presentent une section transversale 
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fractale. La section transversale de chaque pore est constitute suivant un processus 
iteratif (Fig.l). On divise le demi-perimetre d'un cercle de rayon R 0 en T| parties et en 
rernpla9ant chacune de ces rj parties par un demi-cercle ou rainure. A chaque etape k 
du processus, N k nouvelles rainures semi-circulaires de rayon R k et de section totale 
5 Aksont creees. 

La dimension fractale DL de la section transversale a la fin de l'etape k est 
reliee au nombre d'objets Nk genere avec l'echelle donnee I K , par la relation : 

La dimension fractale peut etre deduite d ! une courbe de pression capillaire de 
10 mercure selon la procedure suivante. Du mercure est injecte dans un milieu poreux 
avec une pression d'injection qui croit par paliers. La loi de Laplace permet de 
deduire le volume des pores connaissant le volume de mercure injecte pour une 
pression d'injection donnee et Ton peut construire la courbe de pression capillaire en 
drainage reliant la pression d'injection a la quantite de mercure injectee et la courbe 
15 reliant la proportion du volume total occupe par les pores et la taille des pores. Dans 
le cas ou Ton draine un liquide mouillant hors du milieu poreux tel que de l'eau par 
une injection de gaz, la correlation entre la pression capillaire gaz-eau et la saturation 
de la phase mouillante est donnee par : 

i 

20 Les resultats - experimentaux montrent aisement que les valeurs des 

permeabilites relatives gaz-eau exprimees en fonction des trois saturations, obtenues 
a partir des expressions qu'en donnent les modeles connus et les modes de repartition 
des phases dans la structure des pores, sont loin des valeurs mesurees et done que les 
modeles en cause se revelent par trop simplistes pour representer les interactions 

25 complexes qui se produisent eritre les phases fluides. 
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Par le brevet FR 2 772 483 (US 6 021 662), on connait une methode de 
modelisation pour optimiser de fagon plus rapide et plus realiste les conditions de 
deplacement dans un milieu poreux mouillable par un premier fluide (de l'eau par 
exemple), d'un melange de fluides incluant ce fluide mouillant et au moins un autre 
5 fluide (de l'huile et eventuellement du gaz). Cette methode fait appel a une 
modelisation des pores du milieu poreux par une distribution de capillaires a 
repartition fractale en considerant, dans le cas d'un melange triphasique eau (fluide 
mouillant)-huile-gaz par exemple une stratification des constituants a Tinterieur des 
pores, avec l ? eau au contact des parois, le gaz au centre et l'huile formant une couche 

10 intercalate. Elle comporte une determination experimentale de la courbe de variation 
de la pression capillaire dans les pores en fonction de la saturation relativement aux 
phases liquides, d'ou Ton deduit les valeurs de dimension fractale correspondant a 
une serie de valeurs donnees de la saturation relativement a la phase liquide. Elle 
comporte aussi une modelisation des permeabilites relatives directement sous forme 

15 ^expressions analytiques dependant des differentes valeurs de dimension fractale 
obtenues et en accord avec la repartition stratifiee des differents fluides dans les 
pores. Un simulateur de milieu poreux est utilise a partir de ces permeabilites 
relatives pour determiner les conditions optimales de deplacement des fluides dans le 
milieu poreux. 

20 Le phenomene d'hysteresis fait reference aux variations des proprietes 

petrophysiques (permeabilites relatives, pression capillaire, indice de resistivite, 
etc ...) que Ton observe suivant qu'on precede a la mesure en drainage ou en 
imbibition (ces modes correspondent respectivement & une augmentation et a une 
diminution de saturation de la phase non-mouillante). Ce phenomene doit done etre 

25 pris en consideration pour fournir des valeurs de permeabilites relatives 
, representatives. 

L'etat de la technique concernant les effets d' hysteresis dans les milieux 
diphasiques et triphasiques, est decrit par exemple dans les publications suivantes : 
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- Land C.S. : « Calculation of imbibition relative permeability for two and three- 
phase flow from rock properties », Trans AIME 1968, Vol 243, 149 ; 

- Larsen J. A. , Skauge A. Methodology for numerical simulation with cycle- 
dependent relative permeability", SPEJ, June 1998 ; et 

5 - Carlson F.M. : « Simulation of relative permeability hysteresis to the non wetting 
phase »\ SPE 10157, ATCE, San Antonio Texas, 4-7 Oct 1981. 

La Fig.6 represente typiquement 1' allure des courbes de permeabilites 
diphasiques resultant d'un drainage jusqu'a la saturation irreductible en fluide 

mouillant (M) puis d'une imbibition jusqu'a la saturation residuelle en fluide non- 
10 mouillant (NM). L'hysteresis se manifesto h. deux niveaux. A saturations egales S g , 
on obtient des valeurs numeriques differentes et le point limite final atteint est une 
inconnue qui depend du point de rebroussement S g M a partir duquel on change de 
mode de deplacement. On attribue couramment ce phenomene a la fraction de fluide 
non-mouillant piege. A saturation egale, il n'y a done pas la meme quantite de fluide 
15 mobile ce qui altere les caracteristiques de Pecoulement. 

Pratiquement tous les modeles prenant en compte l'hysteresis font appel £ la 

relation semi-empirique de Land : — — = C L (1) 

S gr S gi 

ou C L represente la constante de Land. Cette relation relie la saturation initiale S gi et 

la saturation residuelle S gr en fluide non mouillant, pour evaluer les saturations en 

20 fluide non mouillant piege et libre. En supposant que cette relation est valable quelle 
que soit la saturation, on 1' applique pour determiner les fractions mobiles 
intermediaires en cours de deplacement. En diphasique, Tassociation de cette relation 
avec des modeles de permeabilites fournit des resultats satisfaisant. 
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En triphasique, l'hysteresis des permeabilites relatives K rg prend une forme 
particuliere. Experimentalement, on observe non seulement une hysteresis de 
deplacement comme c'est le cas en diphasique (importance du sens de variation des 
saturations) mais aussi une hysteresis de cycles puisque les permeabilites sont 
5 dependantes de l'historique des saturations. Sur la Fig.7, les courbes de permeabilites 
relatives K rg correspondant a un premier cycle de drainage et d'imbibition 
(respectivement Dl et II) sont distinctes des courbes correspondahtes (D2, 12) d'un 
deuxieme cycle. 

Le modele de permeabilites relatives de Larsen et al prend en compte ces 
10 deux formes d'hysteresis. Partant d'une approche combinant le modele de Stone en 
parallele avec la formule de Land et la methode ^interpolation de Carlson, approche 
ou seule l'hysteresis de deplacement est prise en compte, il faut le noter, Larsen et al 
ont introduit un facteur de reduction empirique fonction de la saturation en eau qui 
permet d'approcher la reduction de permeabilite du gaz associee k l'hysteresis de 
15 cycles. 

Definition de la methode selon l'invention 

La methode de modelisation selon l'invention permet d f optimiser de fagon 
plus rapide et plus realiste les conditions de deplacement dans un milieu poreux 
mouillable par un premier fluide, de melanges di ou triphasiques incluant ce premier 

20 fluide mouillant et au moins un deuxieme fluide non mouillant. Elle offre de ce fait 
aux operateurs, un outil plus fiable pour evaluer notamment les meilleurs modes de 
deplacement des fluides au sein du milieu poreux. en drainage et en imbibition. Elle 
est basee sur une representation fractale du milieu poreux avec une modelisation des 
pores par une distribution de capillaires a section fractale en considerant une 

25 repartition stratifiee des fluides a Tinterieur des pores, le fluide mouillant se 
repartissant au contact des parois et autour du deuxieme fluide (ou des deux autres en 
triphasique). 
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La methode comporte en combinaison : 

la determination experimental de la courbe de variation de la pression capillaire 
(Pc) dans les pores d'un echantillon de ce milieu poreux en presence d'un fluide 
mouillant et d f au mpins un fluide non mouillant, (par injection de mercure dans 
5 un echantillon place sous vide par exemple) ; 

la determination a partir de cette courbe de pression capillaire, des valeurs de la 
dimension fractale correspondant a une serie de valeurs donnees de la saturation 
relativement aux liquides ; 

- la modelisation des effets d'hysteresis modifiant les saturations mobiles des 
10 fluides effectivement deplaces en fonction du nombre de cycles de drainage et 

d'imbibition subis par 1' echantillon, en faisant intervenir des constantes de 
piegeage ou de depiegeage du fluide non mouillant differentes selon que Ton 
realise une phase de drainage ou d'imbibition ; 

la modelisation des permeabilites relatives directement sous forme ^expressions 
15 analytiques dependant des differentes valeurs de dimension fractale obtenues ; et 

- Introduction des permeabilites relatives dans un simulateur du milieu poreux et 
la determination au moyen de ce simulateur, des conditions optimales de 
deplacement des fluides du melange dans le milieu poreux. 

La methode s'applique par exemple a la determination des deplacements de 
20 melanges de fluides comportant un premier fluide mouillant, un deuxieme fluide non 
mouillant et un gaz, en considerant une repartition stratifiee des fluides a Tinterieur 
des pores, le fluide mouillant se repartissant au contact des parois, le gaz occupant le 
centre des pores et le deuxieme fluide etant reparti sous la forme d'une couche 
annulaire au contact a la fois du gaz et du. premier fluide. 

25 La methode peut s f appliquer notamment a la determination par un simulateur 

de gisement, des caracteristiques optimales de substances additionnees a des 
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bouchons de fluide mouillant injectes dans une formation en alternance avec des 
bouchons de gaz, dans le but de deplacer des hydrocarbures en place, ou bien encore 
celles d ! un fluide injecte dans le sol pour deplacer des substances polluantes. 

La rnodelisation des phenomenes obtenue par la methode presente de grands 
5 avantages. Elle pennet une meilleure adequation avec les resultats obtenus en 
laboratoire car les phenomenes physiques sont mieux pris en compte. Les resultats du 
modele sont meilleurs de ce fait lors d'un changement d'echelle par exemple, pour la 
rnodelisation d'une application sur un champ d'operations. 

Le temps de calcul est reduit compare a celui qui est necessaire quand on 
10 utilise des tables, comme avec les methodes anterieures. La rnodelisation de type 
fractal realise est mieux k meme de traiter les effets d'hysteresis rencontres dans 
Tutilisation des precedes d ! injection de type WAG. 

Les resultats de la methode sont en outre parfaitement integrables dans de 
nombreux simulateurs de reservoir : simulateurs 3D, heterogenes, compositionnels, 
15 etc. 

Sexploitation des resultats par les logiciels ^application est facilitee. II n f est 
plus necessaire de faire des interpolations risquees comme c'est la regie quand on 
opere a partir des valeurs discretes des tableaux de resultats, pour tracer des isoperms 
par exemple. 

20 Presentation des figures 

D ! autres caracteristiques et avantages de la methode selon Tinvention, 
apparaitront h la lecture de la description ci-apres dun exemple non limitatif de 
realisation, en se referant aux dessins annexes ou : 



25 



la Fig.l qui illustre sous forme dun tableau, les liens existant pour un melange 
triphasique entre les valeurs obtenues experimentalement, de permeabilites 
relatives d'un fluide et les saturations pour deux des trois fluides ; 
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- la Fig. 2 montre une representation fractale d'un pore ; 

- la Fig.3 montre de fagon schematique la repartition des phases d'un fluide 
triphasique dans un pore fractal avec Ie fluide mouillant W au contact de la paroi, 
la phase gazeuse G repartie dans la plus grande paitie du volume du pore (rayon 

5 Rg), Thuile O etant une couche entre le fluide mouillant et le gaz ; 

- la Fig.4 montre en fonction de la saturation en eau, la courbe de pression 
capillaire d'un echantillon de gres des Vosges par exemple, dont la pente locale S 
sert a determiner la dimension fractale des pores; et 

- la Fig.5 montre en fonction du temps, les courbes de production obtenues 
10 experimentalement pour du gaz (G) pour de l'eau (W) et pour de l'huile (O), 

comparees aux courbes equivalentes obtenues par simulation en utilisant la 
methode selon l'invention ; 

- les Fig. 6A, 6B montrent les effets de l'hysteresis de deplacement affectant les 
permeabilites relatives du fluide mouillant et du gaz, respectivement K m et ; 

15 - la Fig. 7 montre les effets combines de l'hysteresis de deplacement et de cycles 
observes experimentalement sur la phase gazeuse; 

- la Fig. 8 montre les courbes-cles servant a la modelisation de la constante de 
depiegeage ; 

la Fig. 9 montre le resultat obtenu avec le modele dans le cadre d'une injection de 
20 type WAG ; et 

la Fig 10 montre la validation du modele sur une experience de type WAG pour 
les courbes de recuperation des trois phases. 
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Description detaillee de la methode 

La methode selon l'invention peraiet de determiner les permeabilites relatives 
triphasiques de milieux poreux, en s'appuyant sur un modele de type fractal du milieu 



5 - Kalaydjian, FJ-M et al. : "Three Phase Flow in Water-wet Porous Media: 
Determination of Gas-oil Relative Permeabilities under Various Spreading 
conditions"; 68 th Ann. Tech. Conf. AndExh. of the SPE, Houston, Texas, 1993. 

La methode selon l'invention comporte, comme on l'a vu une modelisation de 
Tecoulement des phases avec repartition des phases h. Tinterieur de la structure 

10 fractale de chaque pore. Dans le cas d'un ecoulement d'eau et d'huile au sein d'une 
structure poreuse mouillable a l'eau, Thuile s'ecoule dans le volume du pore. Dans le 
cas d'un ecoulement triphasique, il y a repartition stratifiee, l'eau qui est le fluide 
mouillant, s'ecoule le long des parois des pores, le gaz circule dans le volume du pore 
et Thuile s'ecoule entre le gaz et l'eau. Les saturations sont calculees comme etant la 

15 surface relative dans une section transversale occupee par chacun des fluides. 

A Tequilibre, toutes les rainures ayant un rayon superieur a Rk qui est donne 
Y 

par la loi de Laplace Pc=2— sont occupees par le gaz, les tubes les plus petits, par 

R k 

les deux autres fluides (eau et huile). Ainsi on exprime la saturation en fluide 
mouillant comme la fraction de l'aire des tubes occupes. 

20 En calculant la fraction de 1'aire des capillaires occupee par Teau pour tous les 

rayons entre R k et R«,., on aboutit a Texpression : 



poreux, en se basant sur une approche decrite par : 
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Y 

et comme Pc=2 — r, la correlation entre la pression capillaire et la saturation de la 
R k 

phase mouillante est donnee par 

p c = s w -^h 

ou S w est la saturation de la phase mouillante. 

5 La representation graphique de cette correlation dans un diagramme log-log 

est une ligne droite partant du point (S w , Pc) correspondant au capillaire le plus grand 
de la structure fractale avec un rayon R 0 . On peut supposer que : 

- le rayon Ro, le premier envahi quand le mercure est injecte (Fig.4), correspond a 
une valeur de la saturation de l'ordre de l/r=10 3 . Chaque segment de la courbe de 

10 pression capillaire est une partie d'une ligne partant de Ro (supposee la meme 

pour tous les differents segments), correspondant a la correlation P c , S w deja 
citee. Chaque ligne a une pente donnee, on peut lui associer une dimension 
lineaire fractale. Les valeurs de la pente varient de -1,5 a -3,3 comme le montre la 
Fig.4, ce qui conduit a des valeurs de la dimension lineaire fractale D L entre 1,3 et 

15 1,7. 

chaque domaine est atteint par le mercure pour des saturations correspondant a 
1'emplacement ou on trouve R 0 sur chaque ligne. 

La saturation des deux liquides quand le gaz est present dans le pore, est 
calculee comme expliquee plus haut pour une phase : 



20 S Lig . = 



El 



2-D L 



en supposant que les deux liquides occupent les tubes dont les rayons sont inferieurs 
ou egaux a R K et le gaz, le centre de chaque pore. La saturation en huile est Taire 
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relative de la section transversale des rainures occupee par lliuile dont le rayon est 
inferieur ou egal a Rk. 

Modelisation de Fhysteresis 

Pour proceder au calcul des permeabilites relatives, il est necessaire de 
5 determiner pour chaque phase la fraction qui circule et par consequent, d'estimer de 
fagon systematique, les saturations correspondant aux fractions stagnantes. Ceci doit 
etre fait pour les deux cas etudies, par exemple le drainage de l'eau et de Thuile par le 
gaz, et rirnbibition en eau. 

L'originalite de la modelisation proposee ici est de prendre en compte 
10 Thysteresis directement a son origine c'est-k-dire au niveau des phenomenes de 
piegeage et depiegeage de la phase non mouillante. Sur la Fig.8, si les courbes 
pouvaient se deduire exactement par translation, cela signifierait qu'une partie du gaz 
est restee piegee au cours du drainage secondaire et ne participe pas a l'ecoulement. 
En realite, ce n'est pas le cas mais le fait qu'a saturation en gaz egale, la penneabilite 
15 soit plus faible au cours de D2 que pendant II signifie que la fraction de gaz piege est 
plus importante au cours du drainage. En d'autres termes, la non reversibilite des 
courbes de permeabilites peut s'expliquer par une non reversibilite (hysteresis) entre 
le phenomene de piegeage et celui de depiegeage. 

On conserve done la formule de Land (equation 1) mais on introduit une 
20 constante de depiegeage valable pendant les phases de drainage, differente de la 
constante de Land valable pendant les phases dimbibition. Tout revient done a 
modeliser revolution de la constante de depiegeage au cours des cycles. 

Les caracteristiques a prendre en compte au niveau de la formulation sont les 
suivantes. On considere que Cp=Cl pendant toutes les phases de piegeage car on se 
25 trouve directement dans les conditions d' application de Land ou une constante 
unique suffit & decrire le phenomene. Pour Sgr faible en fin de drainage, C D est proche 
de C P . On traduit juste ici le fait que le gaz est plus facile a depieger lorsqu'il y en a 
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peu dans la carotte car on peut alors penser qu'il se trouve dans des pores de grosse 
taille facilement reconnectables lors de la reinjection. Quand Sgr augmente, Co 
devient superieur a Q, ce qui traduit un depiegeage moins efficace. 

Pour des saturations en gaz elevees, on retrouve une certaine reversibilite 
5 mais il existe une fraction de gaz piegee inaccessible dans des pores de petite taille. 
On atteint done une courbe de K rg unique basse qui correspond a un cas pseudo 
diphasique ou Thuile n'est plus mobile. 

D'apres cette representation qui repose sur des observations experimentales et 
des hypotheses de travail, la valeur de C D passe par un maximum puisque la 
10 constante de depiegeage est egale a Cl a la fois quand est faible et quand il est 
maximum lorsque Ton decrit la courbe de basse mobilite. L'expression suivante 
permet de concilier tous les aspects precedents : 





( K 1 -K ^ 


E 








rg 1X rgrxun 






(C DM -C L ) + C L 






< s " 2 J 



Cette formulation comporte plusieurs parametres : 
15 Cdm : il prend en compte Tecart entre le piegeage et le depiegeage ; 
Cl : constante de Land ; 
Krgmin : courbe de basse mobilite ; 
E : parametre de calage. 

K 1 ^: valeur de la permeabilite relative au gaz au debut de Timbibition precedente. 

20 K D1 rg : valeur de la permeabilite relative au gaz sur la courbe de premier drainage 
pour la saturation en gaz correspondant a K r rg . 

Quel que soit Tordre et la nature du cycle considere, la relation (2) et la 
relation de Land (1) permettent de determiner les saturations piegees et mobiles par 
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rintermediaire de la formule de Land tout en prenant en compte le caractere 
triphasique de l'hysteresis (deplacement et cycles) sans recours a des facteurs de 
reduction empiriques. 

5 etS g = ^, (4) 

C est egal a C L ou Cd suivant le mode de deplacement ; 
Sgf: saturation en gaz libre ; 
Sgt: saturation en gaz piege. 

Calcul des permeabilites relatives: 

10 A) Permeabilites relatives des liquides 

L'application de la loi de Poiseuille a chaque capillaire du faisceau pour la 
phase qui l'occupe, permet de calculer les permeabilites relatives k l'eau et Hiuile 
(Fig.5). 

Des etudes experimentiales effectuees par ailleurs (Larsen et al) ont montre 
15 qu'il existait une relation permettant de relier la saturation residuelle en huile au 
cours des cycles et la fraction de gaz piege. 

$or = or) Sgt=0'~ a S gt 

- (S 0 r)sgt=o represente la saturation en huile residuelle laissee en place dans le 
milieu avant que du gaz soit piege. 

20 En ne considerant que la fraction en circulation qui contribue a la conductivity 

hydraulique, les permeabilites relatives pour Teau et Thuile s'expriment de la maniere 
suivante: 



(3) 
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R est le facteur de reduction lie au piegeage de la phase non mouillante. 

Dans ces expressions, il est utile de mentionner que : 

5 - la saturation en eau irreductible Swi est supposee stable. 

- La gamme de tailles de capillaires occupee par l'huile mobile est calculee comme 
la difference entre les tailles des capillaires occupes par les deux liquides avec la 
saturation totale en liquide Si=S 0 +S w , et celles des capillaires satures en eau et en 
huile stagnante. 

10 - K ro (2ph.) represente la valeur de la permeabilite relative a l'huile determinee par 
un test d ! imbibition a l'eau et l'huile. Quand les phases eau et huile sont seules 
presentes et comme l'echantillon teste est mouillable a Teau, Thuile s^coulera 
dans la section du pore exactement comme le fait le gaz dans un ecoulement 
triphasique. 

15 B) Permeabilite relative au gaz 

Comme le gaz est une phase non mouillante, il occupe Tespace central du 
pore et il s'etend vers sa peripheric au fur et a mesure que la saturation en gaz 
augmente, sans toutefois arriver au contact de la paroi solide (Fig.6). On considere 
que le gaz circule dans un seul pore dont le rayon R g est donne par la relation : 

20 R g = R 0 +Ri+R 2 + ... + R k , 

la permeabilite au gaz etant donnee alors par : 

a- (Si+^r) 4 (5) 
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avec a = 2=^* > Djl representant la dimension fractale lineaire du milieu poreux et S L 
etant la saturation totale en liquide egale a 1-S g . 

Le modele de peimeabilites relatives permettant de calculer K™, K ro et K rg est 
implante dans un simulateur tel qu ! ATHOS® ou GENESYS©. Cela permet de 
5 proceder au calage d'experiences realisees au laboratoire mais aussi d'optimiser les 
conditions a remplir pour deplacer des fluides petroliers en place dans un gisement, 
que ce soit par injection de gaz soit par injection alternee de bouchons d'eau et de gaz 
(methode dite WAG), en tenant compte des conditions de pression et de temperature 
regnant a la profondeur de production. 

10 Validation 

La methode selon l'invention a ete validee a travers plusieurs types 
d^xperiences : 

On a injecte du gaz dans des milieux poreux contenant de Teau et de Thuile, dans 
des conditions tres variees. On voit sur la Fig. 5 par exemple que Ton obtient un 
15 tres bon accord entre les courbes de production des trois phases (eau, huile, gaz) 

obtenues experimentalement et celles predites par le simulateur de gisement 
alimente par les donnees obtenues en accord avec la methode. 

- On a aussi injecte alternativement du gaz et de l'eau (injection de type WAG). La 
Fig 10 montre que Ton obtient dans ce cas aussi un accord excellent pour les trois 
20 phases tout au long de la duree de Texperience. 
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REVENDICATIONS 

1) Methode de modelisation pour optimiser les conditions de deplacement 
dans un milieu poreux mouillable par un premier fluide, d'un melange de fluides 
triphasique incluant ce premier fluide mouillant et au moins un deuxieme fluide non 
5 mouillant, caracterisee en ce qu'elle comporte en combinaison : 

la determination expeiimentale de la courbe de variation de la pression capillaire 
(P c ) dans les pores d'un echantillon de ce milieu poreux en presence du fluide 
mouillant et du fluide non mouillant ; 

- la modelisation des pores du milieu poreux par une distribution de capillaires a 
10 section fractale en considerant une repartition stratifiee des fluides a Tinterieur 

des pores, le fluide mouillant se repartissant au contact des parois et autour du dit 
au moins un autre fluide ; 

la determination a partir de la elite courbe de pression capillaire (Pc), des valeurs 
de dimension fractale correspondant a une serie de valeurs donnees de la 
15 saturation relativement aux phases liquides ; 

la modelisation des effets d'hysteresis modifiant les saturations mobiles des 
fluides effectivement deplaces dans P echantillon selon le nombre de cycles de 
drainage et d ! imbibition subis par P echantillon, en faisant intervenir des 
constantes de piegeage ou de depiegeage du fluide non mouillant, differentes 
20 selon que Pon est dans un phase de drainage ou une phase d'imbibition ; 

la modelisation des permeabilites relatives directement sous forme ^expressions 
analytiques dependant des differentes valeurs de dimension fractale obtenues ; et 

Introduction des permeabilites relatives dans un simulateur du milieu poreux et 
la determination au moyen de ce simulateur, des conditions optimales de 
25 deplacement du melange de fluides dans le mili eu poreux. 
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2) Methode selon la revendication 1, caracterisee en ce Ton modelise les pores 
du milieu poreux par une distribution de capillaires a repartition fractale en 
considerant une repartition stratifiee des fluides a l'interieur des pores, le fluide 
mouillant se repartissant au contact des parois, le gaz occupant le centre des pores et 

5 le deuxieme fluide etant reparti sous la forme d'un film annulaire au contact a la fois 
du gaz et du premier fluide. 

3) Methode selon la revendication 1 ou 2, caracterisee en ce que Ton met en 
oeuvre le simulateur de gisement pour determiner les caracteristiques optimales de 
substances additionnees a des bouchons de fluide mouillant injectes dans une 

10 formation en alternance avec des bouchons de gaz dans le but de deplacer des 
hydrocarbures en place. 

4) Methode selon Tune des revendications precedentes, caracterisee en ce 
quelle comporte la mise en ceuvre d'un simulateur de gisement pour determiner les 
caracteristiques optimales d'un fluide injecte dans le sol pour drainer des substances 

15 polluantes. 
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